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Abstract. We give a rigorous and self-contained survey of the abelian sand- 
pile model and rotor-router model on finite directed graphs, highlighting the 
connections between them. We present several intriguing open problems. 



1. Introduction 

The abelian sandpile and rotor-router models were discovered several times by 
researchers in different communities operating independently. The abelian sandpile 
model was invented by Dhar |Dha90j as a test-bed for the concept of self-organized 
criticality introduced in [BTW87] . Related ideas were explored earlier by Engcl 
Eng75, Eng76 in the form of a pedagogical tool (the "probabilistic abacus"), by 



Spencer Spe87| pp. 32-35], and by Lorenzini |Lor89[ ILor91j in connection with 



arithmetic geometry. The rotor-router model was first introduced by Priezzhcv 
et al. [PDDK96 (under the name "Eulerian walkers model") in connection with 
self-organized criticality. It was rediscovered several times: by Rabani, Sinclair and 
Wanka [RSW98J as an approach to load-balancing in multiprocessor systems, by 
Propp |Pro01| as a way to derandomize models such as internal diffusion-limited 
aggregation (IDLA) |DF91[|LBG92) . and by Dumitriu, Tetali, and Winkler as part 
of their analysis of a graph-based game [DTW03 ] . Articles on the chip-firing game 
in the mathematical literature include jBig'99; Big97 BLS91, BL92 . Those on the 
rotor-router model include |Lev02| lLP05l IHP08| ILP07a| ILP07bj . Below we briefly 
describe the two models, deferring the formal definitions to later sections. 
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The abelian sandpile model on a directed graph G, also called the chip-firing 
game, starts with a collection of chips at each vertex of G. If a vertex v has at 
least as many chips as outgoing edges, it can fire, sending one chip along each 
outgoing edge to a neighboring vertex. After firing a sequence of vertices in turn, 
the process stops when each vertex with positive out-degree has fewer chips than 
out-going edges. The order of firings does not affect the final configuration, a fact 
we shall discuss in more detail in Section [2] 

To define the rotor-router model on a directed graph G, for each vertex of G, 
fix a cyclic ordering of the outgoing edges. To each vertex v we associate a rotor 
p(v) chosen from among the outgoing edges from v. A chip performs a walk on G 
according to the rotor-router rule: if the chip is at v, we first increment the rotor 
p(v) to its successor e = (v, w) in the cyclic ordering of outgoing edges from v, and 
then route the chip along e to w. If the chip ever reaches a sink, i.e. a vertex of 
G with no outgoing edges, the chip will stop there; otherwise, the chip continues 
walking forever. 

A common generalization of the rotor-router and chip-firing models, the 
height arrow model, was proposed in |PDDK96] and studied in [DR04J . 

We develop the basic theory of the abelian sandpile model in section [2] and 
define the main algebraic object associated with it, the sandpile group of G |Dha90j 
(also called the "critical group" by some authors, e.g. Big99, WagOO ). Further- 
more, we establish the basic results about recurrent chip configurations, which 
play an important role in the theory. In Section [3] we define a notion of recurrent 
configurations for the rotor-router model on directed graphs and give a character- 
ization for them in terms of oriented spanning trees of G. The sandpile group acts 
naturally on recurrent rotor configurations, and this action is both transitive and 
free. We deduce appealing proofs of two basic results of algebraic graph theory, 
namely the Matrix- Tree Theorem |Sta99|, 5.6.8] and the enumeration of Eulerian 
tours in terms of oriented spanning trees )Sta99| Cor. 5.6.7]. We also derive a fam- 
ily of bijections between the recurrent chip configurations of G and the recurrent 
rotor configurations of G. Such bijections have been constructed before, for ex- 
ample in BW97j; however, our presentation differs significantly from the previous 
ones. Section [4] establishes stronger results for both models on Eulerian digraphs 
and undirected graphs. In Section [5] we present an alternative view of the rotor- 
router model in terms of "cycle-popping," borrowing an idea from Wilson's work 
on loop-erased random walk; see [PW98] . We conclude in Section [6] by presenting 
some open questions. 



2. Chip-Firing 

In a finite directed graph (digraph) G = (V, E), a directed edge e € E points from 
the vertex tail(e) to the vertex head(e). We allow self-loops (head(e) = tail(e)) 
as well as multiple edges (head(e) = head(e') and tail(e) = tail(e')) in G. The 
out-degree outdeg(w) of a vertex v (also denoted by d v ) is the number of edges e 
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with tail(e) = v, and the in-degree indeg(u) of v is the number of edges e with 
hcad(e) = v. A vertex is a sink if its out-degree is zero. A global sink is a sink s 
such that from every other vertex there is a directed path leading to s. Note that 
if there is a global sink, then it is the unique sink. 

If G has the same number of edges from v to w as from w to v for all vertices 
v w then we call G bidirected. In particular, a bidirected graph is obtained by 
replacing each edge of an undirected graph with a pair of directed edges, one in 
each direction. 

Label the vertices of G as vi, v^, ■ ■ ■ , v n . The adjacency matrix A of G is the 
n x n matrix whose (i,j)-entry is the number of edges from Vi to Vj, which we 
denote by a VuVj or ay . The (graph) Laplacian of G is the nxn matrix A = D — A, 
where D is the diagonal matrix whose (i, £)-entry is the out-degree of Vi. That is, 



Note that the entries in each row of A sum to zero. If the vertex m is a sink, then 
the i th row of A is zero. 

A chip configuration a on G, also called a sandpile on G, is a vector of non- 
negative integers indexed by the non-sink vertices of G, where <j(v) represents 
the number of chips at vertex v. A chip configuration a is stable if o~{v) < d v 
for every non-sink vertex v. We call a vertex v active in a if v is not a sink and 
o~(v) > d v . An active vertex v can fire, resulting in a new chip configuration a' 
obtained by moving one chip along each of the d v edges emanating from v; that 
is, <j'{w) = a(w) + a vw for all w ^ v and o'(v) = <r(v) — d v + a vv . We call the 
configuration a' a successor of a. 

By performing a sequence of firings, we may eventually arrive at a stable chip 
configuration, or we might continue firing forever, as the following examples show. 

Example 2.1. Consider the complete directed graph on three vertices (without self- 
loops). Then placing three chips at a vertex gives a configuration that stabilizes 
in one move, while placing four chips at a vertex gives a configuration that never 
stabilizes (see Figured]). 

It might appear that the choice of the order in which we fire vertices could 
affect the long-term behavior of the system; however, this is not the ColSC , MS the 
following lemma shows (and Figure [2] illustrates). 

Lemma 2.2 ( |Dha90] . |DF91j ). Let G be any digraph, let oq, a±, . . . , a n be a sequence 
of chip configurations on G, each of which is a successor of the one before, and let 
<7q, a[, . . . , a' m be another such sequence with a' Q = <7q. 

1. If a n is stable, then m < n, and moreover, no vertex fires more times in 
(7q , . . . , a' m than in o~o, . . . , o~ n . 

2. If o~ n and a m are both stable, then m = n, a n — o~' n , and each vertex fires the 
same number of times in both histories. 




di — au fori=j. 
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Figure 1. Some chip configurations eventually stabilize, while 
others never stabilize. 




Figure 2. Commutation of the chip-firing operations. 



Proof. Part 2 is an immediate corollary of part 1, which we now prove. If part 1 
fails, then consider a counterexample with m + n minimal. Let Vi be the vertex that 
fires when o~i-i becomes o~i, and v[ be the vertex that fires when o~' i _ 1 becomes a[. 
Vertex v[ must be fired at some stage (say the i th ) in the sequence of configura- 
tions (To, . . • , £>-„, since a n is stable; say Vi — v[. Then Vi,v\, 112, ■ ■ ■ , Ui+i, ■ ■ ■ ,v n 
is a permissible firing sequence that turns do i n t° °Vi with the same number of fir- 
ings at each site as the unpermuted sequence. The firing sequences v\, . . . , £>i-i, f j+i , 
. . . , v n and v' 2 , v' 3 , . . . , v' m then constitute a smaller counterexample to the lemma 
(with initial configuration a[), contradicting minimality. □ 

Definition 2.3. Starting from a configuration a, Lemma |2~21 shows that there is at 
most one stable configuration that can be reached by a finite sequence of firings 
(and that if such a configuration exists, then no infinite sequence of firings is 
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possible). If such a stable configuration exists we denote it er° and call it the 
stabilization of a. 

Thus far, the presence or absence of sinks was irrelevant for our claims. For 
the rest of this section, we assume that the digraph G has a global sink s. 

Lemma 2.4. If digraph G has a global sink, then every chip configuration on G 
stabilizes. 

Proof. Let N be the number of chips in the configuration. Given a vertex v of G, 
let vq, Ui, . . . , iV-ij Vr be a directed path from vq = v to v r = s. Every time v r —i 
fires, it sends a chip to the sink which remains there forever. Thus w r -i can fire 
at most ./V times. Every time u r -2 fires, it sends a chip to tv_i, and d Vr _ 1 such 
chips will cause v r -\ to fire once, so 1V-2 fires at most d Vr _ 1 N times. Iterating 
backward along the path, we see that v fires at most d Vl ■ ■ ■d Vr _ 1 N times. Thus 
each vertex can fire only finitely many times, so by Lemma 12.21 the configuration 
stabilizes. □ 

We remark that when G is connected and the sink is the only vertex with 
in-degree exceeding its out-degree, the bound one gets from the above argument 
on the total number of firings is far from optimal; see Theorem 14.81 for a better 
bound. 

Define the chip addition operator E v as the map on chip configurations that 
adds a single chip at vertex v and then lets the system stabilize. In symbols, 

E v a={a + l v )° 

where l v is the configuration consisting of a single chip at v. 

Lemma 2.5. On any digraph with a global sink, the chip addition operators com- 
mute. 

Proof. Given a chip configuration a and two vertices v and w, whatever vertices 
are active in a + l v are also active in configuration a' = a + l v + l w . Applying 
to a' a sequence of firings that stabilizes a + l v , we obtain the configuration 
E v a + l u ,. Stabilizing this latter configuration yields E w E v o~. Thus E w E v o~ is a 
stabilization of a'. Interchanging the roles of v and w, the configuration E v E w a is 
also a stabilization of a' . From Lemma 12721 we conclude that E w E v a — E v E w a. □ 

Lemma [2751 is called the abelian property; it justifies Dhar's coinage "abelian 
sandpile model" . From the above proof we also deduce the following. 

Corollary 2.6. Applying a sequence of chip addition operators to a yields the same 
result as adding all the associated chips simultaneously and then stabilizing. 

Let G be a digraph on n vertices with global sink s. The reduced Laplacian A' 
of G is obtained by deleting from the Laplacian matrix A the row and column 
corresponding to the sink. Note that firing a non-sink vertex v transforms a chip 
configuration a into the configuration a — , where A' v is the row of the reduced 
Laplacian corresponding to v. Since we want to view the configurations before and 
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after firing as equivalent, we are led to consider the group quotient Z" 1 /H, where 
H = Z™ _1 A' is the integer row-span of A'. 

Definition 2.7. Let G be a digraph on n vertices with global sink s. The sandpile 
group of G is the group quotient 

S(G) = Z n - 1 /Z n - 1 A'{G). 

The connection between the sandpile group and the dynamics of sandpiles 
on G is made explicit in Corollary 12.161 For the graph in Figure [31 the sandpile 
group is the cyclic group of order 3. The group structure of 5(G) when G is a tree 
is investigated in |Lev07| . 

Lemma 2.8. The order ofS(G) is the determinant of the reduced Laplacian A'(G). 

Proof. The order of S(G) equals the index of the lattice H = Z"- 1 A' in Z"" 1 , 
and, recalling that the volume of a parallelepiped is the determinant of the matrix 
formed from its edge- vectors, we deduce that this in turn equals the determinant 
of A'. □ 

Lemma 2.9. Let G be a digraph with a global sink. Every equivalence class ofTH 1 ^ 1 
modulo A'(G) contains at least one stable chip configuration of G . 

Proof. Let <5 be the configuration given by S(v) = d v for all v, and let 8° be its 
stabilization. Then S°(v) < d v for all v ^ s, so 8 — 8° is a positive vector equivalent 
to the zero configuration. Given any a £ Z n , let m denote the minimum of all 
the coordinates of a together with (so that m < 0). Then the vector 

[3 = a + (-m){8 - 8°) 

is nonnegative and equivalent to a. Hence (3° is a stable chip configuration in the 
same equivalence class as a. □ 

Example 2.10. An equivalence class may contain more than one stable chip config- 
uration. For example, consider the complete directed graph on three vertices, with 
one of the vertices made into a sink by deletion of its two outgoing edges (see Fig- 
ure 12]) ■ It is easy to see that there are two stable configurations in the equivalence 
class of the identity: the configuration in which each of the two non-sink vertices 
has chips and the configuration in which each of the two vertices has 1 chip. It 
might seem natural that, if either of these two configurations is to be preferred 
as a representative of the identity element in the sandpile group, it should be the 
former. However, this instinct is misleading, as we now explain. 

Definition 2.11. A chip configuration a is accessible if from any other chip config- 
uration it is possible to obtain cr by a combination of adding chips and selectively 
firing active vertices. A chip configuration that is both stable and accessible is 
called recurrent. 

Remark 2.12. There are several definitions of "recurrent" that are used in the lit- 
erature. Lemma \'2 . 1 71 below shows that these definitions (including the one above) 
are equivalent for any digraph with a global sink. 
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Figure 3. Two stable chip configurations in the equivalence class 
of the identity. 

We will see shortly (in Lemmas 12.131 and I2.15|) that each equivalence class 
in Z n_1 /Z n_1 A' contains a unique recurrent chip configuration. It is customary 
to represent each element of the sandpile group by its unique recurrent element. 
In Example 12.101 the all-1 configuration is accessible, but the all-0 configuration 
is not. Therefore the all-1 configuration is taken as the canonical representative 
of the identity. (However, in the context of cluster-firing and superstabilization 
as described in Definition 14.31 and Lemma 14.61 the all-0 configuration will be the 
preferred representative.) 

Lemma 2.13. Let G be a digraph with a global sink. Every equivalence class of 
Z n_1 modulo A'(G) contains at least one recurrent chip configuration of G. 

Proof. Given a £ Z™" 1 , let m denote the minimum of all the coordinates of a 
together with 0, so that m < 0. Write <i max for the maximum out-degree of a 
vertex in G. Then a is equivalent to the configuration 

/3 = a+[d max + (-m)](S-S°) 

with 6 as in the proof of Lemma 12.91 Since 8 — 8° has all entries positive, we 
have (3 > [d max + (— m)](S — 8°) > d ma>i (S — 8°) > d max . (Inequalities between 
two vectors or between a vector and a scalar are interpreted componentwise.) In 
particular, (3 is accessible, since any chip configuration can first be stabilized, so 
that each vertex has fewer than d max chips, and then supplemented with extra 
chips to obtain 0. Therefore any configuration obtained from (3 by firing is also 
accessible. In particular the stabilization [3° is thus recurrent and equivalent to 
a. □ 

Next we will show that every equivalence class modulo A' contains at most 
one recurrent configuration, making use of the following lemma. 

Lemma 2.14. Let e = (28) — (28)° , where 8 is given by S(v) = d v as before. If a is 
recurrent, then (a + e)° = a. 

Proof. If a is recurrent then it is accessible, so it can be reached from 8 by adding 
some (non-negative) configuration £ and selectively firing. But since a is also stable 
this implies that (£ + 8)° = a. Consider the configuration 

1 = (( + S) + e = 2S + ( + S-(2S)°. 
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Since e > 0, we may start from 7 and fire a sequence of vertices that stabilizes 
£ + 5, to obtain the configuration a + e. On the other hand, since 6 — (25)° > 
we may start from 7 and fire a sequence of vertices that stabilizes 25, to obtain 
the configuration (26)° + ( + 5 — (2S)° = ( + 5, which in turn stabilizes to a. By 
Lemma [2721 it follows that (a + e)° = a. □ 

Lemma 2.15. Let G be a digraph with a global sink. Every equivalence class of 
7h n ~ x modulo A' '(G) contains at most one recurrent chip configuration of G. 

Proof. Let o~± and 02 be recurrent and equivalent mod A'. Label the non-sink ver- 
tices v-y, . . . , v n -i- Then o\ = a 2 + J2ie.J c i^v wnere the Cj are nonzero constants, 
A^ is the row of the reduced Laplacian A' corresponding to u,, and the index i 
runs over some subset J of the integers 1, . . . , n — 1. Write J = J_ U J + , where 
J- = {i : Ci < 0} and J+ = {£ : Cj > 0}, and let 

a = 0-1 + (-Ci)A{ = * 2 + 53 C * A «- 

Let e denote the everywhere-positive chip configuration defined in Lemma 12.141 
Take k large enough so that a' = a + ke satisfies a'(vi) > \ci\d Vi for all i. Starting 
from a', we may fire each vertex Vi for i G J- a total of — a times, and each 
of the intermediate configurations is a valid chip configuration because all the 
entries are nonnegative. The resulting configuration o~± + ke then stabilizes to o~\ 
by Lemma 12.141 Likewise, starting from a' we may fire each vertex Vi for i G J + a 
total of Ci times to obtain 02 + ke, which stabilizes to 02- By Lemma l2.2l it follows 
that (7i = (T2. □ 

Corollary 2.16. Let G be a digraph with a global sink. The set of all recurrent chip 
configurations on G is an abelian group under the operation (o~,o~') 1— * (a + <r')°, 
and is isomorphic via the inclusion map to the sandpile group S(G). 

Proof. Immediate from Lemmas 12.131 and 12.151 □ 

In view of this isomorphism, we will henceforth use the term "sandpile group" 
to refer to the group of recurrent configurations. 

It is of interest to consider the identity element / of the sandpile group, i.e. 
the unique recurrent configuration equivalent to the all-0 configuration. Here is one 
method to compute /. Let a be the configuration 26 — 2. (Arithmetic combinations 
of vectors and scalars are to interpreted componentwise.) Since a° < 6 — 1 we have 
a — o~° > 6 — 1, so a — cr° is accessible. Since a — a° is equivalent to 0, the identity 
element is given by / = (a — a°)° . 

Figure 3] shows identity elements for the Lx L square grid with "wired bound- 
ary," for several values of L. (To be more precise, the graph G is obtained by re- 
placing each edge of the undirected square grid with a pair of directed edges, and 
adjoining a sink vertex s along with two edges from each of the four corner vertices 
to s and one edge from each of the other boundary vertices to s.) The identity 
element of this graph was studied by Le Borgne and Rossin [B R02] , but most basic 
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properties of this configuration, such as the existence of the large square in the 
center, remain unproved. 




Figure 4. The identity element of the sandpile group of the L x L 
square grid for different values of L, namely L = 128 (upper left), 
198 (upper right), 243 (lower left), and 521 (lower right). The 
color scheme is as follows: orange=0 chips, red=l chip, green=2 
chips, and blue=3 chips. 



Figure [5] shows another example, the identity element for the 100 x 100 di- 
rected torus. (That is, for each vertex € Z/100Z x Z/100Z, there are directed 
edges from (i, j) to (i + 1 mod 100, j) and to (i, j + 1 mod 100), and we make (0, 0) 
(the lower-left vertex) into a sink by deleting its two outgoing edges.) 
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Figure 5. The identity element of the sandpile group of the 
100 x 100 directed torus (left) and the 500 x 500 directed torus 
(right). The color scheme is as follows: white=0 chips, black=l 
chip, and the sink, which is at the lower-left corner, is shown in 
red. 



Figure [5] shows a third example, the identity element for a disk-shaped region 
of Z 2 with wired boundary. Examples of identity elements for graphs formed from 
portions of lattices other than the square grid can be found in LKG90 . 




Figure 6. The identity element of the sandpile group of disk- 
shaped regions of diameter 100 (left), 512 (middle), and 521 
(right). The color scheme is as follows: orange=0 chips, red=l 
chip, green=2 chips, and blue=3 chips. 



Also of interest is the inverse of a recurrent chip configuration a (that is, 
the recurrent chip configuration a such that (a + a) — I). One way to compute 
the inverse is via a = (( — (° — er)°, where C is any chip configuration satisfying 
£ > 36 — 3. (Here ( — — a is accessible, since it has at least d v — 1 chips at each 
vertex v.) 
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Given two chip configurations a and £, we say that a is reachable from £ 
(via excitation-relaxation operations) if there exists a configuration (3 such that 
a = (£+/3)°. Note that this implies that a is stable. A digraph is strongly connected 
if for any two distinct vertices v, w there are directed paths from v to w and from w 
to v. We write G\ s for the graph obtained from G by deleting the vertex s along 
with all edges incident to s. 

Lemma 2.17. Let G be a digraph with a global sink s, and let a be a chip configu- 
ration on G. The following are equivalent. 

(1) a is recurrent; that is, a is reachable from any configuration £. 

(2) If £ is any configuration reachable from a, then a is reachable from £. 

(3) a is reachable from any configuration of the form E v o~, where v is a non-sink 
vertex of G. 

(4) Each strongly connected component of G \ s contains a vertex v such that a 
is reachable from E v o~. 

Proof. Since trivially (1) => (2) => (3) =>■ (4), it suffices to show (4) => (1). 

If (4) holds, there is a chip configuration a such that (a + a)° = a and 
a is nonzero on at least one vertex of each strongly connected component of G. 
There exists a positive integer k such that selective firing from ka results in a chip 
configuration (3 with at least one chip at each vertex. Moreover, (a + f3)° = a. 

Now let £ be any chip configuration. Since (3 has at least one chip at each 
vertex, we have £ < a + i[3 for some integer I. Thus we may add chips to £ and 
then stabilize to obtain the configuration (a +10)° = a. Hence a is recurrent. □ 

We also note that, for a digraph G with a global sink, the sandpile group 
is isomorphic to the additive group of harmonic functions modulo 1 on G that 
vanish on the sink [Sol99j . A function / : V(G) — > [0, 1) is harmonic modulo 1 if 
dvf(v) — X™ a v,wf{w) mod 1 for all vertices v. For a sandpile configuration u, 
the associated harmonic function / is the solution of 

w 

For the graph in Figure^ the three harmonic functions are (f(v i), / (^2)) = (0, 0), 
(f(vx), fM) = (1/3, 2/3), and (f(vi), f{v 2 )) = (2/3, 1/3). 

We conclude this section by pointing out a link between the sandpile group 
and spanning trees. By Lemma 12.81 the order of the sandpile group of G equals 
the determinant of the reduced Laplacian A' of G. By the matrix-tree theorem 
Sta99l 5.6.8], this determinant equals the number of oriented spanning trees of G 
rooted at the sink (that is, acyclic subgraphs of G in which every non-sink vertex 
has out-degree 1). Various bijections have been given for this correspondence; see, 
for example, [BW97 . In Section [3] we will use the rotor-router model to describe a 
particularly natural bijection, and deduce the matrix-tree theorem as a corollary. 
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3. Rotor-Routing 

Chip-firing is a way of routing chips through a directed graph G in such a fashion 
that the chips emitted by any vertex v travel in equal numbers along each of the 
outgoing edges. In order to ensure this equality, however, chips must wait at a 
vertex v until sufficiently many additional chips have arrived to render v active. 
Rotor-routing is an alternative approach to distributing chips through G which 
dispenses with this waiting step. Since we cannot ensure exact equality without 
waiting, we settle for the condition that the chips emitted by any vertex v travel 
in nearly equal numbers along each of the edges emanating from v. We ensure 
that this near-equality holds by using a rotor mechanism to decide where each 
successive chip emitted from a vertex v should be routed. 

Given a directed graph G, fix for each vertex v a cyclic ordering of the edges 
emanating from v. For an edge e with tail v we denote by e + the next edge after e 
in the prescribed cyclic ordering of the edges emanating from v. 

Definition 3.1. A rotor configuration is a function p that assigns to each non-sink 
vertex v of G an edge p(v) emanating from v. If there is a chip at a non-sink 
vertex v of G, routing the chip at v (for one step) consists of updating the rotor 
configuration so that p(v) is replaced with p(v) + , and then moving the chip to 
the head of p(v) + . A single-chip-and-rotor state is a pair consisting of a vertex w 
(which represents the location of the chip) and a rotor configuration p. The rotor- 
router operation is the map that sends a single-chip-and-rotor state (w,p) (where 
w is not a sink) to the state (w + ,p + ) obtained by routing the chip at w for one 
step. (See Figure [7] for examples of the rotor-router operation.) 

As we will see, there is an important link between chip-firing and rotor- 
routing. A hint at this link comes from a straightforward count of configurations. 
Recall that a stable chip configuration is a way of assigning some number of chips 
between and d v — 1 to each non-sink vertex v of G. Thus, the number of stable 
configurations is exactly Y[ v 4, where the product runs over all non-sink vertices. 
This is also the number of rotor configurations on G. Other connections become 
apparent when one explores the appropriate notion of recurrent states for the 
rotor-router model. We will treat two cases separately: digraphs with no sink, and 
digraphs with a global sink (Lemma 13.61 applies to both settings) . 

Definition 3.2. Let G be a sink-free digraph, i.e. one in which each vertex has at 
least one outgoing edge. Starting from the state (w, p), if iterating the rotor-router 
operation eventually leads back to (w, p) we say that (w, p) is recurrent; otherwise, 
it is transient. 

Our first goal is to give a combinatorial characterization of the recurrent 
states, Theorem l3.8l We define a unicycle to be a single-chip-and-rotor state (w, p) 
for which the set of edges {p(v)} contains a unique directed cycle, and w lies on 
this cycle. (Equivalently, p is a connected functional digraph, and w is a vertex on 
the unique cycle in p.) The following lemma shows that the rotor-router operation 
takes unicycles to unicycles. 
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Lemma 3.3. Let G be a sink-free digraph. If (w, p) is a unicycle on G, then (w + , p + ) 
is also a unicycle. 

Proof. Since (w, p) is a unicycle, the set of edges {p(v)} v ^ w = {p + (v)} v ^ w contains 
no directed cycles. The set of edges {p + (v )} forms a subgraph of G in which every 
vertex has out-degree one, so it contains a directed cycle. Since any such cycle 
must contain the edge p + {w) = p(w) + , this cycle is unique, and w + lies on it. □ 

Lemma 3.4. Let G be a sink-free digraph. The rotor-router operation is a permu- 
tation on the set of unicycles of G. 

Proof. Since the set of unicycles is finite, by Lemma 13.31 it is enough to show 
surjectivity. Given a unicycle U — (u>,p), let U~ = (w~ , p~) be the state obtained 
by moving the chip from w to its predecessor w~ in the unique cycle through w, 
and replacing the rotor at w~ with its predecessor in the cyclic ordering of outgoing 
edges from w~ . Then the rotor-router operation applied to U~ yields U. It remains 
to show that U~ is a unicycle; for this it suffices to show that every directed cycle 
in p~ passes through w~ . Suppose that there is a directed cycle of rotors in p~ 
which avoids w~. Since p~ agrees with p except at w~ , this same directed occurs 
within p and avoids w~ , a contradiction since w~ is on p's unique cycle. □ 

Corollary 3.5. Let G be a sink-free digraph. If {w, p) is a unicycle on G, then (w, p) 
is recurrent. 

In Lemma [3.71 below, we show that the converse holds when G is strongly 
connected. We will need the following lemma, which is analogous to Lemma [2.41 
for the abelian sandpile. A vertex w is globally reachable if for each other vertex 
v there is a directed path from v to w. 

Lemma 3.6. Let G be a digraph with a globally reachable vertex w. For any start- 
ing vertex and rotor configuration, iterating the rotor-router operation a suitable 
number of times yields a state in which the chip is at w. 

Proof. Since w is globally reachable, either G is sink-free or w is the unique sink. 
Thus either we can iterate the rotor-router operation indefinitely, or the chip even- 
tually visits w. In the former case, since G is finite, the chip visits some vertex 
v infinitely often. But if ir is a vertex that is visited infinitely often and there is 
an edge from x to y, then y is also visited infinitely often. Inducting along a path 
from v to w, we conclude that the chip eventually visits w. □ 

Lemma 3.7. Let G be a strongly connected digraph. If (w, p) is a recurrent single- 
chip- and-rotor state on G, then it is a unicycle. 

Proof. Since G is strongly connected, every vertex is globally reachable. Hence 
by Lemma 13.61 if we start from any initial state and iterate the rotor-router rule 
sufficiently many times, the chip visits every vertex of G. 

Suppose (w, p) is a recurrent state. Once every vertex has been visited and we 
return to the state (w, p), suppose the rotors at vertices Vi . . . , Vk form a directed 



14 



Holroyd, Levine, Meszaros, Peres, Propp and Wilson 



cycle. If w does not lie on this cycle, then for each i, the last time the chip was 
at Vi it moved to and hence the edge from to was traversed more 

recently than the edge from to Uj. Carrying this argument around the cycle 
leads to a contradiction. Thus, every directed cycle in the rotor configuration must 
pass through w. But now if we start from w and follow the rotors, the first vertex 
we revisit must be w. Hence (w,p) is a unicycle. □ 

Combining Corollary 13.51 and Lemma 13. 7\ we have proved the following. 

Theorem 3.8. Let G be a strongly connected digraph. Then (w,p) is a recurrent 
single- chip- and-rotor state on G if and only if it is a unicycle. 

Next we consider the case when G is a digraph with a global sink. Note that 
we cannot apply the rotor-router operation to states in which the chip is at the 
sink. We call these absorbing states. For any starting state, if we iterate the rotor- 
router operation sufficiently many times, the chip must eventually arrive at the 
sink by Lemma 13.61 

A chip-and-rotor state is a pair r = (a, p) consisting of a chip configuration a 
and rotor configuration p on G. A non-sink vertex is active in r if it has at least 
one chip. If v is active, then firing v results in a new chip-and-rotor state given by 
replacing the rotor p(v) with p{v) + and moving a single chip from v to the head of 
p{v) + (and removing the chip if p(v) + is a sink). We say that t' is a successor of r if 
it is obtained from r by firing an active vertex. We say that r is stable if no vertex 
can fire, i.e., all chips have moved to a sink and disappeared. The rotor-router 
operation has the following abelian property analogous to Lemma 12.21 

Lemma 3.9. Let G be a digraph with a global sink. Let tq, t\, . . . , r„ be a sequence 
of chip-and-rotor states of G, each of which is a successor of the one before. If 
To, r{, . . . , t^j is another such sequence, and t„ is stable, then m < n. If in addition 
T' m is stable, then m — n and t„ = r' n , and for each vertex w, the number of times 
w fires is the same for both histories. 

Proof. Let Vi and v[ be the vertices that are fired in t^_i and t[_ 1 to obtain Ti 
and t[, respectively. We will show that if r„ is stable and the sequences v and 
v' agree in the first i — 1 terms for some i < m, then some permutation of v 
agrees with v' in the first i terms. Since v[ is active in Ti_i = T 'i-ii it must be 
active in Ti, t^i, . . . , until it is fired. Since T n is stable, it follows that Vj = v[ 
for some j > i. Let j be the minimal such index. Starting from to, the vertices 
Vi, V2, ■ • ■ , Vi-i,Vj,V{, Vi+x, . . . , Vj-i, Uj'+i, . . . ,v n can be fired in that order, result- 
ing in the same stable configuration r n . Moreover, this sequence agrees with v' in 
the first i terms. 

By induction, it follows that the sequence v' is an initial subsequence of a 
permutation of v. In particular, m < n. If T' m is also stable, by interchanging the 
roles of r and r , we obtain that v' is a permutation of v. □ 
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Figure 7. The unicycle configurations resulting from the evolu- 
tion of a particular unicycle on the bidirected 3-by-4 rectangular 
grid. By Lemma l4~9l the chip traverses each directed edge exactly 
once before the original unicycle is obtained. Thus the number 
of distinct unicycle configurations equals the number of directed 
edges, in this case 34. From Lemma [4.111 it follows that from any 
given unicycle, after some number of steps, the state will be the 
same but with the cycle's direction reversed. This occurs, for ex- 
ample, with unicycles 1 and 13. 



Given a non-sink vertex v in G, the chip addition operator E v is the map 
on rotor configurations given by adding a chip at vertex v and iterating the rotor- 
router operation until the chip moves to the sink. By Lemma l3.9l and the reasoning 
used in the proof of Lemma 12.51 the operators E v commute. This is the abelian 
property of the rotor-router model. 
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If, rather than running the chips until they reach the sink, each chip is run for 
a fixed number of steps, then the abelian property fails, as the example in Figured] 
illustrates. (The proof of Lemma 13.91 requires that chips be indistinguishable, and 
it is not possible to run each chip for a fixed number of steps without distinguishing 
between them.) Despite the failure of commutativity, this way of routing chips has 
some interesting properties, similar to the bound given in Proposition 13.211 see 
work of Cooper and Spencer [CS06J. 




Figure 8. Failure of the abelian property for rotor-router walk 
stopped after two steps. 

A rotor configuration p on G is acyclic if the rotors do not form any directed 
cycles. If G has a global sink, then p is acyclic if and only if the rotors form an 
oriented spanning tree rooted at the sink. 

Lemma 3.10. Let G be a digraph with a global sink, and let v be a vertex of G. The 
chip addition operator E v is a permutation on the set of acyclic rotor configurations 
on G. 

Proof. We first argue that applying E v to an acyclic rotor configuration yields 
an acyclic rotor configuration: this is proved by induction on the number of rotor- 
routing steps, where the induction hypothesis states that every path leads to either 
the sink or to the chip. 

Since the set of acyclic rotor configurations is finite, it suffices to show sur- 
jectivity. Let p be an acyclic rotor configuration on G, add an edge e from the 
sink s to v to form a sink-free digraph G' , and assign the rotor p(s) — e. Then 
U = (s,p) is a unicycle on G' . Starting from U, we can iterate the inverse of the 
rotor-router operation (which by Lemma l3~4l is well-defined for unicycles) until the 
next time we reach a state U' = (s,p') with the chip at s. Because U' is a unicycle 
with chip at s, deleting the edge e leaves an acyclic rotor configuration p' . Observe 
that running the rotor-router operation from U' to U, upon ignoring the edge e, 
is equivalent to applying E v to p' and obtaining p. □ 

Definition 3.11. We next describe an action of the sandpile group on acyclic rotor 
configurations. Given a chip configuration a and a rotor configuration p on G, write 
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o~(p) for the rotor configuration obtained by adding u{v) chips at each vertex v 
and routing them all to the sink. By Lemma 13.91 the order of these routings is 
immaterial. Thus we may write cr(p) as 

*(p)= f n k (v) )p, 

\vev(G) J 

where the product symbol represents composition of the operators. 
It is trivial that <T2{<Ji{p)) — (<ti + 02) (p). 

Since acyclic rotor configurations on G can be identified with oriented span- 
ning trees rooted at the sink, Lemma 13.101 implies that every chip configuration a 
acts as a permutation on the set of oriented spanning trees of G rooted at the sink. 

Lemma 3.12. Let G be a digraph with a global sink, and let p be an acyclic rotor 
configuration on G. If the chip configurations g\ and o~2 o,re equivalent modulo the 
reduced Laplacian A' of G, then o~\(p) = o~2(p). 

Proof. If cr(v) > d v , and we route (for one step) d v of the chips at v, then the rotor 
at v makes one full turn and one chip is sent along each outgoing edge from v. By 
Lemma T3.91 it follows that if a' is a successor to a (that is, o' is obtained from a 
by firing some active vertex v), then <r(p) — cr'(p) for any rotor configuration p. 
Inducting, we obtain o~(fl) = cr°(p) for any rotor configuration p. 

In particular, if / is the recurrent chip configuration that represents the iden- 
tity element of the sandpile group, we have 

I(I(p)) = (I + I)(p) = (I + I)°(p)=I(p) 

for any rotor configuration p. By Lemma l3.10l thc map p t— ► J(p) is a permutation 
on the set of acyclic rotor configurations, so it must be the identity permutation. 
Now if <7i , (72 are equivalent modulo A', then (<ti + 1) and {a 2 + 1)° are recurrent 
and equivalent modulo A', hence equal by Lemma \2. 151 Since 

o-i(p) = o-i(I(p)) = (at + I)(p) = & + I)°(p), i = 1,2, 

we conclude that <j\ (p) = 02 (p) ■ □ 

It follows from Lemma 13.121 that the sandpile group of G acts on the set of 
oriented spanning trees of G rooted at the sink. Our next lemma shows that this 
action is transitive. 

Lemma 3.13. Let G be a digraph with a global sink. For any two acyclic rotor 
configurations p and p' on G, there exists a chip configuration a on G such that 
o{p)=p>. 

Proof. For a non-sink vertex v, let a(y) be the number of edges e such that p(v) < 
e < p'(v) in the cyclic ordering of outgoing edges from v. Starting with rotor 
configuration p, and with a(v) chips at each vertex v, allow each chip to take 
just one step. The resulting rotor configuration is p'; let (3 be the resulting chip 
configuration, so that a(p) = P(p'), and let 7 be the inverse of the corresponding 
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element (j3 + 1) of the sandpile group. By Lemma T3. 121 and the fact that (3 + 7 is 
equivalent to modulo A', we have 



Next we define recurrent rotor configurations on a digraph with a global sink, 
and show they are in bijection with oriented spanning trees. 

Definition 3.14. Let G be a digraph with a global sink. Given rotor configurations p 
and p' on G, we say that p is reachable from p 1 if there is a chip configuration a 
such that a(p') = p. We say that p is recurrent if it is reachable from any other 
configuration p' . 

Note that in contrast to Definition 13. 2\ the location of the chip plays no role 
in the notion of recurrent states on a digraph with global sink. 

Lemma 3.15. Let G be a digraph with a global sink. A rotor configuration p on G 
is recurrent if and only if it is acyclic. 

Proof. By Lemma l3.10[ any configuration reachable from an acyclic configuration 
must be acyclic, so recurrent implies acyclic. Conversely, if p is acyclic and p' is 
any rotor configuration, the configuration l(p') (where 1 denotes the configuration 
with one chip at each vertex) is acyclic, since the rotor at each vertex points along 
the edge by which a chip last exited. By Lemma 13 . 131 there is a chip configuration 
a such that o~(l(p')) — p, so p is reachable from p' and hence recurrent. □ 

Just as for the sandpile model, there are several equivalent definitions of 
recurrence for the rotor-router model. 

Lemma 3.16. Let G be a digraph with a global sink s, and let p be a rotor configu- 
ration on G. The following are equivalent. 

(1) p is acyclic. 

(2) p is recurrent; that is, p is reachable from any rotor configuration p' . 

(3) // p' is reachable from p, then p is reachable from p' . 

(4) p is reachable from any rotor configuration of the form E v p, where v is a 
vertex of G. 

(5) Each strongly connected component of G \ s contains a vertex v such that p 
is reachable from E v p. 

Proof. By Lemma GUI we have (1) (2), and trivially (2) (3) (4) (5). 

If property (5) holds, let C±, . . . ,C(, be the strongly connected components 
of G \ s, and for each i, let Vi £ Ct be such that p is reachable from E Vi p. Choose 
an integer k large enough so that if we start k chips at any Vi and route them to 
the sink, every vertex in Ci is visited at least once. Let 



(a + 7 )(p) = ([3 + 7 )(p') = p'. 
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Then in p' ', the rotor at each vertex points along the edge by which a chip last 
exited, so p' is acyclic. Since p is reachable from p' , by Lemma f3. 101 it follows that 
p is acyclic. Thus (5) (1), completing the proof. □ 

Next we show that the action of the sandpile group on the set of oriented 
spanning trees of G is free. 

Lemma 3.17. Let G be a digraph with a global sink, and let <j\ and <ji be recurrent 
chip configurations on G. If there is an acyclic rotor configuration p of G such that 
eri(p) = o- 2 (p), then o x = a 2 - 

Proof. Let a = 0-1+02 (recall that 02 is the inverse of 02.) By Lemma 13.121 
o-(p) = ^(o-i(p)) = CT2(cr 2 (p)) = ((72 +f?2~)(p) = P- Since cr(p) = p, after adding 
a to p, for each vertex v, the rotor at v makes some integer number c v of full 
rotations. Each full rotation results in d v chips leaving v, one along each outgoing 
edge. Hence a = J2v Cv^, which is in the row span of the reduced Laplacian, so 
a is equivalent to modulo A'(G), and hence a\ — 02 by Lemma T2.15I □ 

Corollary 3.18 (Matrix Tree Theorem). Let G be a digraph and v a vertex of G. 
The number of oriented spanning trees of G rooted at v is equal to the determinant 
of the reduced Laplacian A'(G) obtained by deleting from A(G) the row and column 
corresponding to v. 

Proof. Without loss of generality we may assume the graph is loopless, since loops 
affect neither the graph Laplacian nor the number of spanning trees. 

If v is not globally reachable, then there are no spanning trees rooted at v, 
and there is a set of vertices S not containing v, such that there are no edges in 
G from S to S c . The rows of A'(G) corresponding to vertices in S sum to zero, so 
A'(G) has determinant zero. 

If v is globally reachable, delete all outgoing edges from v to obtain a di- 
graph G' with global sink v. Note that G and G' have the same reduced Laplacian, 
and the same set of oriented spanning trees rooted at v. 

Fix an oriented spanning tree p of G'. The mapping a 1— » u{p) from S{G') 
to the set of oriented spanning trees of G' is a surjection by Lemma 13.131 and is 
one-to-one by Lemma [3371 and by Lemma l2~8l \S{G')\ = det A'(G')- □ 

Given a digraph G with global sink, define its rotor-router group as the 
subgroup of permutations of oriented spanning trees of G generated by the chip 
addition operators E v . 

Lemma 3.19. LL07] The rotor-router group for a digraph G with a global sink is 
isomorphic to the sandpile group S(G). 

Proof. The action of the sandpile group on oriented spanning trees is a homo- 
morphism from the sandpile group S(G) onto the rotor-router group. For any 
two distinct sandpile group elements o~\ and cr 25 f° r an y oriented spanning tree 
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p (there is at least one), by Lemma 13.171 we have <J\{p) ^ c 2 (p), so the associ- 
ated rotor-router group elements are distinct, i.e., the group homomorphism is an 



Since the number of recurrent chip configurations of G equals the number 
of oriented spanning trees, it is natural to ask for a bijection. Although there is 
no truly "natural" bijection, since in general there is no canonical spanning tree 
to correspond to the identity configuration, we can use the rotor-router model to 
define a family of bijections. Fix any oriented spanning tree p rooted at the sink, 
and associate it with the identity configuration i\ For any other oriented spanning 
tree p', by Lemma 13.131 there exists a £ S(G) with a(p) = p'\ moreover, a is 
unique by Lemma 13.171 Associate a with p' . Since this defines a surjective map 
from recurrent configurations to oriented spanning trees, it must be a bijection. 

Remark 3.20. A variant of the rotor-router rule relaxes the cyclic ordering of 
edges emanating from a vertex, and merely requires one to choose some edge 
emanating from the current location of the chip as the new rotor-setting and move 
the chip along this edge. This is the branching operation introduced by Propp 
and studied by Athanasiadis Ath97 . Alternatively, one can put a probability 
distribution on the edges emanating from each vertex, and stipulate that the new 
edge is to be chosen at random. This gives the tree- walk introduced by Anantharam 
and Tsoucas |AT89j in their proof of the Markov chain tree theorem of Leighton 
and Rivest |LR86j . 

We conclude this section with the following result of Holroyd and Propp 
HP08 , which illustrates another area of application of the rotor-router model. 

Proposition 3.21. Let G = (V,E) be a digraph and let Y C Z be sets of vertices. 
Assume that from each vertex there is a directed path to Z. Let a be a chip config- 
uration on G. If the chips perform independent simple random walks on G stopped 
on first hitting Z, let H(o~,Y) be the expected number of chips that stop in Y. If 
the chips perform rotor-router walks starting at rotor configuration p and stopped 
on first hitting Z , let H p (a, Y) be the number of chips that stop in Y . Then 



Here e~ is the edge preceding e in the cyclic ordering of edges emanating from u. 
Since h is a harmonic function on G, the sum of h{u) — h(v) over all edges e = (u, v) 
emanating from u is zero, so the formula wt(e) = h{u) — h{v) + wt(e~) remains 



isomorphism. 



□ 



\H p (a,Y) -H(a,Y)\ < ^ |ft(head of e) - ft (tail of e)\ 



(1) 



edges e 
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valid even when e = p{u). We assign weight wt(p(i>)) to a rotor configuration p, 
and weight h(v) to a chip located at v. By construction, the sum of rotor and 
chip weights in any configuration is invariant under the operation of rotating the 
rotor at a chip and then routing the chip. Initially, the sum of all chip weights is 
H(a,Y). After all chips have stopped, the sum of the chip weights is H p (a,Y). 
Their difference is thus at most the change in rotor weights, which is bounded 
above by the sum in (JTJ). □ 

Similar bounds hold even for some infinite directed graphs in which the right 
side of {!]) is not finite. Thus rotor-routing can give estimates of hitting probabil- 
ities with very small error. See |HP08j for more details. 



4. Eulerian Graphs 

A digraph G = (E, V) is Eulerian if it is strongly connected, and for each vertex 
v e V the in-degree and the out-degree of v are equal. We call G an Eulerian 
digraph with sink if it is obtained from an Eulerian digraph by deleting all the 
outgoing edges from one vertex; equivalently, G has a sink and every other vertex 
has out-degree that is at least as large as its in-degree. An Eulerian tour of a 
digraph G is a cycle in G that uses each edge exactly once. Such a tour exists 
if and only if G is Eulerian. Note that for any connected undirected graph, the 
corresponding bidirected graph is Eulerian. In this section we show some results 
that do not hold for general digraphs, but are true for Eulerian ones. We first treat 
the sandpile model, and then the rotor-router model. 

Lemma 4.1 (Burning algorithm [Dha90'J). Let G be an Eulerian digraph with sink. 
A chip configuration a is recurrent if and only if (a + (3)° — a, where 

f3(v) — outdeg(w) — indeg(w) > 0. 

If a is recurrent, each vertex fires exactly once during the stabilization ofo~ + [3. 

Proof. By the "(4) =>• (1)" part of Lemma l2.17[ if (a + f3)° = a, then a is recurrent. 
Conversely, suppose a is recurrent. Label the non-sink vertices v±, . . . , v n —±. Since 

n-l 

= J2 A 'i> ( 2 ) 

the configurations a and (<7 + J3)° are both recurrent and equivalent modulo A'. 
By Lemma 12.151 it follows that they are equal. 

Let Ci be the number of times vertex fires during the stabilization of a + (3. 

Then 

n-l 
i=l 

The rows of A' are linearly independent, so from ^ we deduce = 1 for alH. □ 
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Informally, the burning algorithm can be described as follows: to determine 
whether a is recurrent, first "fire the sink" to obtain the configuration a + (3. Then 
a is recurrent if and only if every non-sink vertex fires in the stabilization of a + (3. 
In the non-Eulerian case, there is a generalization of the burning algorithm known 
as the script algorithm, due to Eugene Speer |Spe93| . 

Let H be an induced subgraph of G not containing the sink. We say that H 
is ample for a chip configuration a on G if there is a vertex v of H that has at 
least as many chips as the in-degree of v in H . 

Lemma 4.2. Let G be an Eulerian digraph with sink s. A stable chip configuration a 
on G is recurrent if and only if every nonempty induced subgraph of G \ s is ample 
for a. 

Proof. If a is recurrent, there is a chip configuration a such that (5 + a)° = a, 
where S(v) = d v . Each vertex of G fires at least once in the process of stabilizing 
5 + a. Given a nonempty induced subgraph H of G, let v be the vertex of H which 
first finishes firing. After v finishes firing, it must receive at least as many chips 
from its neighbors as its in-degree in H, so o~(v) is at least the in-degree of v in H . 
Thus H is ample for a. 

Conversely, suppose that every nonempty induced subgraph of G is ample 
for a. Let (3 be the chip configuration defined in Lemma 14. II Starting from a + (3, 
fire as many vertices as possible under the condition that each vertex be allowed 
to fire only once. Let H be the induced subgraph on the set of vertices that do not 
fire. Since each vertex v of H is unable to fire even after receiving one chip from 
each incoming edge whose other endpoint lies outside H, we have 

a(v) + d v — indeg H (v) < d v — 1. 

Thus H is not ample and consequently must be empty. So every vertex fires once, 
after which we obtain the configuration a + f3 — Y^,7=i ^'i = "■ Hence (a + f3)° = a, 
which implies a is recurrent by Lemma 12.171 □ 

Next we define a variant of chip- firing called cluster- firing (see Figure |9]), and 
we use Lemma l4~2l to characterize the stable states for cluster-firing. This gives rise 
to a notion of "superstable states" which are in some sense dual to the recurrent 
states. 

Definition 4.3. Let G be a digraph with a global sink. Let a be a chip configuration 
on G, and let A be a nonempty subset of the non-sink vertices of G. The cluster- 
firing of A yields the configuration 

If a' is nonncgative we say that the cluster A is allowed to fire. We say that a is 
superstable if no cluster is allowed to fire. 

Note that a cluster A may be allowed to fire even if no subset of A is allowed 
to fire. For example, in the first configuration in Figure [9l a cluster of two vertices 
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is allowed to fire even though the configuration is stable, so no single vertex is 
allowed to fire. 




stable stable superstable 



Figure 9. A sequence of cluster-firings resulting in a superstable 
chip configuration. The bottom vertex is the sink. The clusters 
that fire are first the two neighbors of the sink, next the top 
vertex, and finally all three non-sink vertices. 



One could consider an even more general operation, "multicluster-firing," 
in which different vertices can be fired different numbers of times, so long as at 
the end of the firings, each vertex has a nonnegative number of chips. However, 
this further-generalized firing operation does not yield anything new for Eulerian 
digraphs, since any multicluster-firing can be expressed as a sequence of cluster- 
firings: Let to denote the maximal number of times that a vertex fires in the 
multicluster-firing, and Cj denote the set of vertices that fire at least j times in 
the multicluster-firing. Since the digraph is Eulerian, C m may be cluster-fired, and 
so by induction the sets C m , C m -u ■ • • jCi can be cluster-fired in that order. 

Denote by 8 the chip configuration 5(v) — d v in which each vertex has as 
many chips as outgoing edges, and by 1 the configuration with a single chip at 
each vertex. 

Theorem 4.4. Let G be an Eulerian digraph with sink. A chip configuration a on G 
is superstable if and only if 5 — 1 — a is recurrent. 

Proof. A cluster A is allowed to fire if and only if for each vertex v £ A we have 

a(v) — d v + indeg A (w) > 0. 

This is equivalent to d v — 1 — o{v) < indeg^(w), i.e., the induced subgraph on A 
is not ample for 5 — 1 — a. By Lemma 1431 the proof is complete. □ 

By Lemmas 12.131 and 12.151 every equivalence class modulo A' contains a 
unique recurrent configuration, so we obtain the following. 

Corollary 4.5. Let G be an Eulerian digraph with sink. Every eguivalence class 
modulo A'(G) contains a unique superstable configuration. 

As a consequence, we obtain that the cluster-firing model on Eulerian di- 
graphs is abelian; this was proved by Paoletti [Fao07J in the bidirected case. 
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Figure 10. Stable sandpilc (on left) and superstable sandpile (on 
right) of 100,000 chips, obtained by placing 100,000 chips at the 
origin of the integer lattice Z 2 and (super)stabilizing. The color 
scheme is as follows: white=0 chips, red=l chip, green=2 chips, 
and blue=3 chips. 



Corollary 4.6. Let G be an Eulerian digraph with sink. Let <jq, a±,...,a n be a 

sequence of chip configurations on G, each of which is obtained from the one before 
by a cluster- firing, with a n superstable. Then any sequence of cluster-firings that 
starts from <7o and ends in a superstable configuration ends in a„ . 

We call the configuration a n in Corollary 14.61 the superstabilization of Co- 

The following result provides a way to compute the superstabilization. 

Proposition 4.7. Let a be a chip configuration on an Eulerian digraph with sink. 
The superstabilization of a is given by 

a* = 6 - 1 - (S - 1 - a° + I)° 

where I is the identity element of the sandpile group. 

Proof. Since the configuration £ = (5 — 1 — a° + 1)° is reachable from the identity 
element, it is recurrent, hence a* = S — 1 — £ is superstable by Theorem 14. 41 Since 
a and a* are equivalent modulo A', it follows from Corollary 14.51 that a* is the 
superstabilization of a. □ 

Our final result concerning the sandpile model on Eulerian digraphs is a 
theorem of Van den Heuvel [vdHOlj ; see also |Tar88| . We give a shorter and more 
direct proof than that presented in vdHOl]. By an bidirected graph with sink s, 
we will mean the digraph obtained from an undirected graph by first replacing 
each edge by a pair of directed edges in opposite directions, and then deleting all 
outgoing edges from s. The effective resistance between two vertices of G is an 
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important quantity in electrical network theory; see, e.g., [DS84] . In particular, the 
quantity i? max appearing in the proposition below is always bounded above by the 
diameter of G, but for many graphs it is substantially smaller than the diameter. 

Proposition 4.8. Let G be a bidirected graph with sink, and let a be a chip config- 
uration on G. The total number of chip moves needed to stabilize a is at most 

2m\a\ R max 

where m is the number of edges, \o~\ is the total number of chips, and i? max is the 
maximum effective resistance between any vertex of G and the sink. 

Note that firing a vertex v consists of d v chip moves. 

Proof. Let c = (To, 01, . . . , ak = o° be a sequence of chip configurations with <7j_|_i 
obtained from cr.; by firing a single active vertex xi . Define the weight of Oi to be 

wt(o-j) = 2J cri(x) wt(x), 

X 

where 

wt(x) = E X T S 

is the expected time for a simple random walk started at x to hit the sink. By 
conditioning on the first step X\ of the walk, we compute 

Awt(a) = E x (E Xl T s - T.) = — 1, 

so firing the vertex xi decreases the total weight by d Xi . Thus 

Wt((7j) - Wt((7j + l) = d Xi . (3) 

By CRR+97 , the function wt is bounded by 2rni? max . Since the final weight 
wt(er°) is nonnegative, summing ([3]) over i we obtain that the total number of chip 
moves TV needed to stabilize a is at most 

fc-i 

N = ^d Xi = wt(cr) - wt(cr°) < 2m \a\ i? max . □ 

i=0 

Next we present results about the rotor-router model specific to the Eulerian 
case. An example of the next lemma is illustrated in Figure [7] 

Lemma 4.9. Let G be an Eulerian digraph with m edges. Let U — (w, p) be a 

unicycle in G. If we iterate the rotor-router operation m times starting from U , 
the chip traverses an Eulerian tour of G, each rotor makes one full turn, and the 
state of the system returns to U . 

Proof. Iterate the rotor-router operation starting from U until some rotor makes 
more than a full turn. Let it be the rotor at vertex v. During this process, v must 
emit the chip more than d v times. Hence if v ^ w, then v must also receive the 
chip more than d v times. Since G is Eulerian, this means that some neighboring 
vertex u must send the chip to v more than once. However, when the chip goes 
from u to v for the second time, the rotor at u has executed more than a full turn, 
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contradicting our choice of v. Thus when the rotor at w has made a full turn, the 
rotors at the other sites have made at most a full turn. 

We can now repeat this argument starting from the configuration obtained 
after the rotor at w has made a full turn. In this way, the future history of the 
system is divided up into segments, each of length at most m, where the chip is 
at w at the start of each segment. It follows that over the course of the future 
history of the system, the chip is at w at least d w /m of the time. 

Since G is strongly connected, we may apply this same argument to every 
state in the future history of the system, with every vertex of G playing the role 
of w. As the system evolves, the chip is at v at least d v /m of the time. Since 
d v /m = 1, the chip is at v exactly d v /m of the time. Hence, as the rotor at w 
executes a full turn, the rotors at the other sites also execute a full turn. Since 
every rotor makes a full turn, every edge is traversed exactly once, so the chip 
traverses an Eulerian tour. □ 

We can use Lemma 14.91 to give a bijective proof of a classical result in enu- 
merative combinatorics relating the number of Eulerian tours of an Eulerian graph 
G to the number of oriented spanning trees of G (see, e.g., |Sta991 Cor. 5.6.7]). 

Corollary 4.10. Let G = (V, E) be an Eulerian digraph. Fix an edge e E E and let 
tail(e) = w. Let T(G,v) denote the number of oriented spanning trees in G rooted 
at w, and let e(G, e) be the number of Eulerian tours in G starting with the edge e. 
Then 

e(G,e)=T(G,w) JJ (d„ - 1)!. 

Proof. There are n«ev(^ — !)■ ways to fix cyclic orderings of the outgoing edges 
from each vertex. There are T(G, w) ways to choose a unicycle U = (w, p) with the 
chip at w and the rotor p(w) = e~, where e~ is the edge preceding e in the cyclic 
ordering of outgoing edges from w. Given these data, we obtain from Lemma 14.91 
an Eulerian tour of G starting with the edge e, namely the path traversed by the 
rotor-router walk in m steps. 

To show that this correspondence is bijective, given an Eulerian tour starting 
with the edge e, cyclically order the outgoing edges from each vertex v in the order 
they appear in the tour. Let p(w) = e~ and for v ^ w let p(v) be the outgoing 
edge from v that occurs last in the tour. Then U — (w,p) is a unicycle. □ 

The following result was first announced in [PPS98J. in the case of rotor- 
router walk on a square lattice. 

Corollary 4.11. Let G be a bidirected planar graph with the outgoing edges at each 
vertex ordered clockwise. Let (w,p) be a unicycle on G with the cycle C oriented 
clockwise. After the rotor-router operation is iterated some number of times, each 
rotor internal to C has performed a full rotation, each rotor external to C has not 
moved, and each rotor on C has performed a partial rotation so that C is now 
oriented counter-clockwise. 
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Proof. Let G' be the graph obtained from G by deleting all vertices and edges 
external to C. Note that G', like G, is Eulerian. Let p~ be the rotor configuration 
on G 1 obtained from p by "regressing" each rotor whose vertex lies on C; that is, 
if v lies on C and p(v) is the edge e, let p(v) = e~ , the edge immediately preceding 
e in the cyclic ordering of the edges emanating from v. Consider the unicycle 
U~ = (w,p~) on G' and note that it has C oriented counter-clockwise. Starting 
from U~ and applying the rotor-router operation #C times, the chip will traverse 
the cycle C , resulting in the state U = (w, p\g')- By Lemma POl further iteration 
of the rotor-router operation on G' returns the system to the state U~. Since the 
outgoing edges at each vertex are ordered clockwise, it is straightforward to see 
that applying the rotor-router rule to U on G' and to (w, p) on G results in the 
same evolution up until the time that state U~ on G' is reached. □ 

Lemma 4.12. Let G be an Eulerian digraph, and let G v be the Eulerian digraph 
with sink obtained by deleting the outgoing edges from vertex v. Then the abelian 
sandpile groups S(G V ) and S(G W ) corresponding to different choices of sink are 
isomorphic. 

Proof. Recall that the sandpile group S(G V ) is isomorphic to Z n_1 /Z" _1 A'(G); we 
argue that for Eulerian digraphs G it is also isomorphic to Z n /Z n A(G). Vectors in 
Z" -1 are isomorphic to vectors in Z™ whose coordinates sum to 0, and modding out 
a vector in Z™ -1 by a row of the reduced Laplacian A' corresponds to modding out 
the corresponding vector in Z™ by the corresponding row of the full Laplacian A. 
For Eulerian digraphs G, the last row of the full Laplacian A is the negative of the 
sum of the remaining rows, so modding out by this extra row has no effect. □ 

We mention one other result that applies to undirected planar graphs, due 
originally to Berman |Ber86| Prop. 4.1]; see also [CROP] . 

Theorem 4.13. If G and G* are dual undirected planar graphs, then the sandpile 
groups of G and G* are isomorphic. (By Lemma \4-.1S\ the locations of the sink are 
irrelevant.) 

5. Stacks and Cycle- Popping 

Let G be a digraph with a global sink. In this section we describe a more general 
way to define rotor-router walk on G, using arbitrary stacks of rotors at each vertex 
in place of periodic rotor sequences. To each non-sink vertex v of G we assign a 
bi-infinite stack p(v) = (pk(v))k£Z of outgoing edges from v. To pop the stack, we 
shift it to obtain (pk+i(v))k£Z- To reverse pop the stack, we shift it in the other 
direction to obtain (pk— i(«))fcez- The rotor-router walk can be defined in terms of 
stacks as follows: if the chip is at vertex v, pop the stack p(v), and then move the 
chip along the edge pi(v). We recover the ordinary rotor-router model in the case 
when each stack p{v) is a periodic sequence of period d v in which each outgoing 
edge from v appears once in each period. 
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The collection of stacks p — (p(v)), where v ranges over the non-sink vertices 
of G, is called a stack configuration on G. We say that p is infinitive if for each 
edge e = (v,w), and each positive integer K, there exist stack elements 

Pk(v) = Pk'( v ) = e 

with k > K and k! < —K. This condition guarantees that rotor-router walk 
eventually reaches the sink. 

Given a stack configuration p, the stack elements pa{v) define a rotor con- 
figuration on G. We say that p is acyclic if po contains no directed cycles. If 
C = {vi, . . . ,v m } is a directed cycle in po, define Cp to be the stack configuration 
obtained by reverse popping each of the stacks p{vi); we call this reverse popping 
the cycle C. (If C is not a directed cycle in pq, set Cp = p.) 

Theorem 5.1. [PW98 Let G be a digraph with a global sink, and let p° be an 

infinitive stack configuration on G. There exist finitely many cycles C\, . . . ,C m 
such that the stack configuration 

p = C m - ■ -Cip° 

is acyclic. Moreover, if C[, . . . ,C' n is any sequence of cycles such that the stack 
configuration p' = C' n ■ ■ ■ C[p° is acyclic, then p' = p. 

If v is a non-sink vertex of G, the chip addition operator E v applied to the 
infinitive stack configuration p is the stack configuration p' obtained by adding a 
chip at v and performing rotor-router walk until the chip reaches the sink. The 
next lemma shows that these operators commute with cycle-popping. 

Lemma 5.2. Let G be a digraph with a global sink, let p be an infinitive stack 
configuration on G, and let C be a directed cycle in G. Then 

E v {Cp)=C{E vP ). 

Proof. Write p' = E v p. Let v — Vq,v\, . . . ,v n — s be the path taken by a chip 
performing rotor-router walk from v to the sink starting with stack configuration p. 
If this path is disjoint from C, then the chip performs the same walk starting with 
stack configuration Cp, and the cycle C is present in p' if and only if it is present 
in po, so the proof is complete. 

Otherwise, choose k minimal and £ maximal with v^,vi £ C. The rotor p' (vi) 
points to a vertex not in C, so the cycle C is not present in p' Q . Thus we must 
show E v (Cp) = p' . With stack configuration Cp, the chip will first travel the path 
vq, . . . ,Vk, next traverse the cycle C, and finally continue along the remainder of 
the path Vk, ■ ■ ■ , v n . Thus the stack at each vertex w G C is popped one more time 
in going from Cp to E v (Cp) than in going from p to E v p; the stack at each vertex 
w C is popped the same number of times in both cases. □ 

The next lemma uses cycle-popping to give a constructive proof of the in- 
jectivity of the chip addition operators E v on acyclic stack configurations. In the 
case of periodic rotor stacks, we gave a non-constructive proof in Lemma 13.101 
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Lemma 5.3. Let G be a digraph with a global sink. Given an acyclic infinitive stack 
configuration p onG and a non-sink vertex v, there exists an acyclic infinitive stack 
configuration p' such that E v p' = p. 

Proof. Let p° be the stack configuration obtained from p by reverse popping the 
stack at each of the vertices on the unique path in po from v to the sink. A rotor- 
router walk started at v with stack configuration p° will travel directly along this 
path to the sink, so E v p° = p. If p° is acyclic, the proof is complete. Otherwise, 
by Theorem 1 5 . 1 1 there are cycles C\, .. . ,C m such that p 1 = C m ■ ■ - Cip° is acyclic. 
By Lemma l5.2[ we have 

E v p' = Cm ■ ■ ■ Ci {E v p°) = C m ■ ■ ■ C\p = p 

where in the last equality we have used that p is acyclic. □ 

Note that the proof shows the following: if p. p' are acyclic infinitive stack 
configurations and E v p' = p, then the unique path in po from v to the sink is 
the loop-erasure of the path taken by rotor-router walk started at v with initial 
configuration p' . 

6. Conjectures and Open Problems 

In this section we discuss some natural questions about chip-firing and rotor- 
routing that remain unanswered. 

Fey-den Boer and Redig [FR07 consider aggregation in the sandpile model 
on Z rf . In their setup, the underlying graph for the chip- firing game is the infinite 
undirected (i-dimensional cubic lattice Z d . Start with each site containing h < 
2d — 2 chips. Here h may be even be taken negative, corresponding to starting 
with a "hole" of depth H = —h at each lattice site; that is, each site absorbs the 
first H chips it receives, and thereafter fires every time it receives an additional 
four chips. Now add n chips to the origin and stabilize. Denote by S Ui h the set of 
sites in 1 d which fired in the process of stabilization. 

Theorem 6.1 ( FR07]). Let C(r) denote the cube of side length 2r + 1 centered at 
the origin in 7L d . For each n there exists an integer r n such that S n .2-2d = C(r n ). 

In two dimensions, Theorem 16.11 states that <SVi,-2 is a square. Simulations 
indicate that for general H > —2, the limiting shape of S n ,H m <^ 2 may be a 
polygon with AH +12 sides. 

Question 6.2. Ln Z 2 , is the limiting shape of S n ^n as in oo a regular (AH + 
\2)-gon? Simulations indicate a regular (AH + \2)-gon with some "rounding" at 
the corners; it remains unclear if the rounded portions of the boundary become 
negligible in the limit. Even if the limiting shape is not a polygon, it would still be 
very interesting to establish the weaker statement that it has the dihedral symmetry 

D 4H -+12- 
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depths H = -2 (left), H = -1 (center), and H = 0. 



The square, octagon and dodecagon corresponding to the cases H — —2, —1,0 
are illustrated in Figure ITT1 Regarding Question 16. 21 we should note that even the 
existence of a limiting shape for S n> u has not been proved in the case H > —2. On 
the other hand, as H — > oo the limiting shape is a ball in all dimensions, as shown 
by Fey and Redig [FR07] and strengthened in |LP07bj . In the theorem below, ujd 
denotes the volume of the unit ball in M d , and B r denotes the discrete ball 

B r = {x e Z d | x\ + . . . + x\ < r 2 }. 

Theorem 6.3. [LP07b| . Fix an integer H > 2 — 2d. Let S n .H be the set of sites in 
1 d that fire in the process of stabilizing n particles at the origin, if every lattice 
site begins with a hole of depth H . Write n — oJd,r d . Then 

Bcir — c-2 C S n ,H 

where 

ci = {2d-l + Hy 1 ' d 
and C2 is a constant depending only on d. Moreover if H > 1 — d, then for any 
e > we have 

S n ,H C 5 c ^ r _|_ c ' 

where 

c[ = {d-e + H)- 1 ^ 
and c' 2 is independent of n but may depend on d, H and e. 

In particular, note that the ratio c\jc\ \ 1 as H f oo. 

For many classes of graphs, the identity element of the sandpile group has 
remarkable properties that are not well understood. Let /„ be the identity element 
of the n x n grid graph G n with wired boundary; the states /„ for four different 
values of n are pictured in Figure [4] Comparing the pictures of /„ for different 
values of n, one is struck by their extreme similarity. In particular, we conjecture 
that as n — > oo the pictures converge in the following sense to a limiting picture 
on the unit square [0, 1] x [0, 1]. 
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Conjecture 6.4. Let a n be a sequence of integers such that a n f oo and — J. 0. For 
x e [0,1] x [0,1] let 

fn{x) = ~2 In ^- 
| \y—nx\ | <a 7l 

There is a sequence a n and a function f : [0, 1] x [0, 1] — > R>o which is locally 
constant almost everywhere, such that f„ — > / at all continuity points of f . 

Most intriguing is the apparent fractal structure in the conjectural /. Recent 
progress has been made toward understanding the fractal structure of the identity 
element of a certain orientation of the square grid; see |CPS07j . 

By Lemma [4.91 the recurrent orbits of the rotor-router operation on an Euler- 
ian digraph are extremely short: although the number of unicycles is typically ex- 
ponential in the number of vertices, the orbits are all of size equal to the number 
of edges. One would expect that such short orbits are not the norm for general 
digraphs. 

Question 6.5. Does there exist an infinite family of non-Eulerian strongly connected 
digraphs G n , such that for each n, all the unicycles of G n lie in a single orbit of 
the rotor-router operation? 

Another question stemming from Lemma [4.9l is the following. Fix two edges eo 
and ei of a digraph G. Starting from a unicycle on G, record a each time the chip 
traverses the edge eo, and record a 1 each time it traverses e\. If G is Eulerian, then 
Lemma T4.9I implies that the resulting sequence will simply alternate 0, 1, 0, 1, . . .. 
For a general digraph, the sequence is periodic, since the initial unicycle must 
recur; what can be said about the period? 

Lastly, the articles jPDDK96j and PPS98] contain several conjectures that 
are supported by both credible heuristics and computer experiments, but that have 
not been rigorously proved. For instance, it appears that, with random initial rotor 
orientations, the set of sites visited by a rotor-router walk of length n in the plane 
typically has diameter on the order of n 1 / 3 [PDDK96] (compare this with the 
corresponding growth rate for random walk in the plane, which is rt 1 / 2 ). 
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